Field theoretic approach to the counting problem of Hamiltonian cycles of graphs 
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A field theoretic representation for the number of Hamiltonian cycles of graphs is studied. By 
integrating out quadratic fluctuations around the saddle point, one obtains an estimate for the num- 
ber which reflects characteristics of graphs well. The accuracy of the estimate is verified by applying 
it to 2d square lattices with various boundary conditions. This is the first example of extracting 
meaningful information from the quadratic approximation to the field theory representation. 
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I. INTRODUCTION 

Let G = (V, E) be a graph with the set of vertices 
V = {rj} and of the edges E = {e^}. A Hamiltonian 
cycle of a graph is a closed path which visits each of the 
vertices once and only once. I denote the number of all 
the Hamiltonian cycles of a graph G by H(G): 



H 



H(G) 



E 



1. 



(1) 



Hamiltonian cycle on G 



See Fig.|l] for examples. 

Hamiltonian cycles have often been used to model col- 
lapsed polymer globules jlj]. The quantity H(G) corre- 
sponds to the entropy of a polymer system in G in a col- 
lapsed but disordered phase. One can model even more 
realistic polymers by introducing a weight that depends 
on the shape of cycles in (|I|) . The polymer melting prob- 
lem is studied by taking into account the bending energy 
with a weight which depends on the number of turns on 
the cycle 0. In ref. q], the protein folding problem is 
studied by incorporating the Van der Waals potential as 
well as the bending energy in the model ([!]). 

For homogeneous graphs (lattices) with N vertices, one 
expects that H(G) behaves like 



H(G) -> C(G)iV T 
where w is defined by 



(TV^oo), 



logw = lim — log H(G). 

N—>oc ly 



(2) 



(3) 



The quantity u> is supposed to be a universal bulk quan- 
tity whereas C(G) and 7 depend on the detail of graphs 
e.g. boundary conditions [Q. 

A field theory representation for ([!]) for arbitrary 
graphs is introduced in ref. [gj and has been used to study 
the extended models [^,^| . For homogeneous graphs with 
the number of vertices N and the coordination number 
q, the saddle point approximation to the representation 
yields 



iV 



(4) 



or oj ~ q/e. This approximation has been proved to 
be good in many examples ||. For a square lattice, 
q/e = 4/e = 1.4715 ... is quite near to the exact value 
u> ~ 1.473 estimated by the direct enumeration |tJ and 
other methods @,|,§. 

In the saddle point approximation, however, graphs 
with identical (g, N) are not distinguished. Indeed there 
is a variety of graphs which has (g, N) in common. Given 
a graph with a pair (q,N), it is often possible to change 
its boundary conditions to modify the 'topology' and the 
'moduli' of the graph keeping the pair. Fig. |l| shows 
examples. Moreover, there are graphs which have identi- 
cal (q, N) but have distinct local structures, e.g. the 2d 
triangular lattice and the 3d cubic lattice. 

In this article, I go beyond the saddle point approxi- 
mation. I work out the quadratic approximation and find 
an estimate for H(G) whose G-dependence is not merely 
through (g, N) but is more sensible. To demonstrate the 
validity of the approximation, I examine 2d square lat- 
tices with a variety of boundary conditions and aspect 
ratios. It is also examined if the estimate to ~ q/e is 
improved. 



II. FIELD THEORETIC REPRESENTATION 

The problem of calculating H(G) is mapped into one 
in a lattice field theory on G H. This is done by in- 
troducing an 0(n) lattice field <j>{r) — (</>i(r), . . . , (pnir)) 
living on V = {rj} with an action 

S[fo)] = \ E &(r)(iA- 1 + e) rr /^(r'). (5) 

The N x N matrix A is the adjacency matrix |l(J of the 
graph G: 

. _ j 1 if r, r' £ V is connected by an e G E 

< II otherwise ' • 



1 



and an infinitesimal parameter e > is introduced for 
convergence. 

The integer H(G) is related to a 2iV-point function by 



H(G) = lim lim - 

n-i-0 n 



Zi 
Z 



Z = D<t> e 



-SWr)] 



Z x = I D<f> e 



-S[<Kr)] 



n 

rev 



</, 2 (r) 



(7) 
(8) 
(9) 



where D$ = Y\ reVA<j < n d<pj(r). Eq. (Q) holds for arbi- 
trary G with N > 2 hecause each terms in the diagram- 
matic expansion corresponds to a Hamiltonian cycle. 



III. APPROXIMATION 



First I evaluate (^) by the saddle point method. I con- 
centrate on Zx since Z — > 1 as n — ► 0. When a graph G 



is homogeneous, there are mean field saddle points which 
are degenerate on 

{j\P(r) = -2iq e }~ (10) 



where 



1 _ 1 

<7e ' Q 



(11) 



This yields the estimate 

Then I consider fluctuations around the saddle point. 
It can easily be seen that there are zero modes corre- 
sponding to the global 0(n)/0(n— 1) symmetry and sub- 
lattice symmetries || . Therefore I am led to introduce a 
gauge fix condition 



5>i(r) = (2<j< 



rev 



By the standard Fadeev-Popov method, I have 



(12) 



Z 1= D<j>e 



•*E r .r' 6 v.i<y£» (r)(iA- 1 +*) PP /*j(r') 



r( 



n u^ n 



2 (r) 



2 ; 2<j<n \reV 



rev 



(13) 



The factor | £ r <M'")r _X 7 I "V r (f ) is the Fadeev- 
Popov determinant multiplied by the volume of the gauge 
orbit. Actually, this is nothing but the jacobian for the 
n-dimensional radial coordinate for constant modes. 

Expanding the field cj> around the saddle point as 

<j)j (r) = y/-2iqJ N , n 5jx + i>j (r), (14) 



where 



one obtains 



IN,: 



n — 1 
2N ' 



(15) 



Zx = 



qjN,. 



NI N , 



2~N r 



7T 2 



x /" ^ e -iE f y 6 v*'W^(^M x e^E^.K^W'l^t")*^) x e-^« 7 ) x [] S (Z>i(»o) , (16) 

2<j<n VrSV / 

where ^4 L (n) and A T (n) are the inverse propagators of the longitudinal mode (j — 1) and the transverse modes 
(2 < j < n), respectively: 



A h rr ,{n) = »(A ^rr' + ^AT^rr' + (1 + I N ] n )^ rr , + —1(1- Z"^), 



iC(n) = t(A ^r-r' - -I N ]jrr> + (1 - /^ n )e<5 rr /. 



(17) 
(18) 



One neglects the interaction terms Vint(V') anc l performs the gaussian integrations to obtain 

' NI N „ n T , JV-i "I 



(2iV) 



r(f) 



det'i A T ( 



det^ A L ( 



x2, 



(19) 



2 



where e is sent to +0. The last factor two corresponds 
to the residual symmetry Y^ r ( Pi( r ) ~ J2r- 4>i( r )- The 
prime on det means the omission of the eigenvalue for 
the constant mode. 

One notices that the signature of the real part (the last 
term) of ( |l~8|) changes at n = 1. Thus one has to assume 
n > 1 to derive (19). Then in eq.(|l9|) one takes the limit 
n — > to obtain the final result 



IV. SQUARE LATTICES 



H(G) 



AT, 



e 2 ( J Ar,o) 



NI N , 



it det^(A T (0)A) 
N det^(A L (0)A) ' 



(20) 



Note that det A det' A T (0) = gdet'(A T (0)A). Eq. © 
is nothing but the estimate I desired to have. The sad- 
dle point result (q/e) N is corrected by the ratio of de- 
terminants which contains information of details of the 
structure of G. 

In ref . j^] , it is claimed that the quadratic correction to 
(0) vanishes. In the present analysis, the Fadeev-Popov 
method is worked out to find y/w/N missing in ref. ||. 
As shown in Sec. IV, the ratio of determinants is not equal 



to unity and contributes to 7 and C(G) non-trivially. 
Moreover, inclusion of i = \/— 1 and e > in the action 
(0) in the present analysis enables one to discuss the limit 
of application of the quadratic approximation in Sec. IV. 



To see how eq.(gOj) works, I study concrete examples 
P(L 1 ,L 2 ) and SP{L ly L 2 ) shown in Fig.|l|. Both are two- 
dimensional square lattices with the edge lengths L\ and 
L 2 . The difference between P and SP lies in boundary 
conditions. For P(L\, L 2 ), the periodic boundary condi- 
tion is imposed for both two directions. For SP{L\,L 2 ), 
that across the edge L 2 is replace by the skew-periodic 
one. They are good examples to test eq.(|o|). One can 
switch the boundary condition or vary the aspect ratio 
L 2 /Li to make the graph globally distinct while keeping 
(q,N) = (4, L\ x L 2 ). The saddle point approximation 
cannot see the difference among them but eq.(^o|) has a 
chance to distinguish them. 

Graphs P(Li,L 2 ) and SP{L\,L 2 ) can be viewed as 
discrete tori with different moduli parameters. Thus it is 
interesting also in its own right to determine the asymp- 
totic behaviors of H(P(L 1 ,L 2 )) and H{SP{L 1 ,L 2 )) in 
the limit L\, L 2 — ► 00. 

I analytically evaluate ( p?i| ) for P{L\,L 2 ) and 
SP(L±, L 2 ). In the momentum representation, the de- 
terminants become (a = for P and a — 1 for SP) 



det'(A T (n)A) = J] 
det(A L (n)A) = 2 x Jj' 



0<n,<L,--l 



<7lj<Lj-l 



— cos ki + 0— ) + cos k 2 

1 ( f, ki 1 

— cos fci + a— ) + cos fc 2 

2iAr,n V V L\ 



(21) 
(22) 



The constant mode k\ = k 2 = is excluded in Y[ ■ Hereafter, the indices rij and kj should be related by kj — 2nn,j/Lj 
The formula 



n 



0<m<i-l 



X — 2x COS 



2m7r 



+ 1 



„2i 



2x L cos(L6) + 1 



(23) 



enables one to obtain 

det , (A T (0)A) 
det(A L (0)A) 

where 

u±(k) -- 



2LiL 2 I Ll L 2y0 Y\ 



u-(k 2 



n 2 =0 



u+(k 2 



y(k) L i 
y(k) L ^ 



y(k)~ Ll ±2, for P, 
y(k)~ Ll - 2cosfc for SP, 



y(k) = 2I LlL2 . 



cos k 



(24) 



(25) 



L 2 ,0 



cos/c) 2 - 1. (26) 



Now I specialize to the case where both L\ and L 2 are 
odd. In the limit N = L\ x L 2 — > 00 with R — L 2 /L\ 
fixed, it is allowed to approximate ]T by exp(J dk 2 log) 
after separating quickly oscillating factors. In this limit, 
assuming < R < 1 for periodic case, I obtain 



HiP^L^J- 



it sin( ^)iexp( 1 ^ r - 
cosh 2 (^-^)^ 



H(SP(ii,L 2 ))~(|)%F 



sin(^)2 sinM^- ^)5 



sinh 2 (^-^)* 



(28) 



One sees the quadratic approximation predicts the ex- 
plicit form of C(G) and the fact 7 = 1 in (g). It is 
remarkable that the correction is independent of N in 
the limit. Namely, it implies that the estimates for u> are 
unchanged. 

FigsB and m show the plot of p7]), ( p8| ) normalized by 



(4/ 



\L{1 



as functions of L 2 jL\ (solid line). The finite 



(27) 



L\L 2 results (^0) for odd Lj such that 3 < Lj < 29 are 
also shown (solid circle) . For the purpose of comparison, 
I plot the exact numbers of Hamiltonian cycles (box) and 
estimates of them by simulations (solid box) for Lj odd 

The exact numbers are determined by direct enumer- 
ation by a program in C. Table 1 shows the result. 
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This program is useful only for small graphs because the 
time needed for calculation increases exponentially with 
L\ x L 2 - The simulation is based on the biased Monte 
Carlo method with a code in C. The time needed grows 
again exponentially with L\ x L 2 but with a smaller expo- 
nent. I have used direct enumeration for L\ x L 2 < 39 and 
Monte Carlo simulation for L\ x L 2 ^ 90. The evalua- 
tion of (^p| ) takes time proportional to L\ x L 2 . Therefore 
eq.(|2^) has an advantage even if it is approximate. 

Figs ^ and || suggest that the quadratic approximation 
to the field theory is reliable. The field theory succeeds 
in predicting that the correction depends almost only on 
L2/L1 and that there is a qualitative difference between 
P and SP. The L 2 /Li-dependence of the correction gen- 
erally agrees with the exact result though there is a slight 
deviation in the small L2/L1 region. 

There is a definite discrepancy at L2/L1 = 7r 2 /2 for 
the skew-periodic case. The quadratic approximation 
( p8| ) diverges at L 2 /Li = tt 2 /2 as seen in Fig.| while 
the exact result takes finite values and is simply increas- 
ing there. Actually, one can argue that L2/L1 > ir 2 /2 
is out of range of application of (|20|). As I mentioned, I 
defined the limit n — > as the continuation from n > 1 
where the gaussian integral converges. Let us look at the 
evolution of the spectrum of (A L (n)A) on the way from 
n > down to n = 0. The eigenvalue for kj = ( 1 — 1 / L j ) ir 
mode hits zero at some < n c < 1 if 

7T TT 

2/i!L 2 ,o < cos — — + cos — . (29) 
L\L2 L2 

This condition is equivalent to L 2 /Li > ir 2 /2 in the limit 
Li x L2 — > 00. So one suddenly has a zero mode at n = n c 
and no symmetry is responsible for it. This suggests that 
the quadratic approximation breaks down there and that 
(p8|) is not reliable for such values of L2/L1. A careful 
analysis shows that the periodic case is free from such 
spurious zero modes. 



V. SUMMARY AND DISCUSSIONS 

I have found an approximate formula for the number 
of Hamiltonian cycles of graphs by the quadratic approx- 
imation to the field theoretic representation. It has been 
tested for 2d square lattice with a variety of boundary 
conditions. I have obtained a multiplicative correction 
which is independent of the size N = #V. Namely, the 
dependence on the boundary conditions is obtained and 
is shown to be in a good agreement with the true behav- 
ior. 

A natural extension to the present problem is the 
counting of the number of closed self-avoiding walks that 
visit fN vertices of a graph, where < / < 1 is fixed. In 
that case, the form (||) can still be assumed but (7(G), 7 
and u> now depend on /. There is a field theoretic rep- 
resentation for u)(f) called lattice cluster theory p2| . It 
is not identical with (Q) even for / = 1. It is interesting 



to compare thequadratic approximation to lattice clus- 
ter theory fL^JTif with the present result. In the former 
one, estimate for uj(f) is improved over the saddle point 
approximation for / < 1. For / = 1, the estimate for 
u> = is unchanged in accord with the present case. 
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FIG. 1. Examples of graphs and Hamiltonian cycles on 
them. One does not distinguish the base points and the di- 
rections of a cycle, (a) H(G) = 1 for this graph, (b) P(3,4), 
the 2d square lattice with the periodic boundary condition. A 
Hamiltonian cycle is drawn in the thick line. For this graph, 
(q,N) = (4, 12). (c) SP(3,4), the 2d square lattice with the 
skew-periodic boundary in the horizontal direction. It is lo- 
cally isomorphic to P(3, 4) but the boundary condition makes 
it distinct globally. 



FIG. 2. The number of Hamiltonian cycles for P(L\,Li) 
for odd Lj as a function of L2/ L\. Normalized by the saddle 
point result (4/e) Llia . Plotted are the quadratic approxima- 
tion to the field theory ( solid square: 3 < Li, L2 < 29, dashed 
line: the limit L\ x L2 — * co ), exact results obtained by enu- 
meration (box, Li X L 2 < 39), and estimates by weighted 
Monte Carlo simulations (solid box, L\ x L2 < 90). 
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enumeration 
Monte Carlo 

field theory (saddle point) j 
field theory (fluctuation, L h L 2 <30) < 
field theory (fluctuation, L^-x^) /s 



■" A 



0.1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

arctan(L 2 /Li) (radian) 

FIG. 3. The log plot of the number of Hamiltonian cycles 
for SP(L\, L2) for odd Lj as a function of arctan(L2/ii)- 
Normalized by the saddle point result (4/e) ili2 . Plotted are 
the same as in Fig.H. 
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□ L-'. enumeration 
A Monte Carlo 

■ ■" /• ' field theory(saddle point) 
/field theory(fluctuation, i,,L 2 <30) 
.....-*" field theory(fluctuation, LixL2->°°) 
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0.4 0.6 

Li/L, 



0.1 



TABLE I. Exact number of Hamiltonian cycles for L\ x Li 
2d square lattice with periodic (P) and skew-periodic (SP) 
boundary conditions. Determined by the direct enumeration. 



Lx 


L 2 


H{P{L X ,L 2 )) 


H(SP(L U L2)) 


H{SP(L 2 ,Li)) 


3 


3 


48 


55 


55 


5 


3 


390 


397 


866 


5 


5 


23580 


29001 


29001 


7 


3 


2982 


2989 


13021 


7 


5 


1045940 


1108006 


1820582 


9 


3 


23646 


23653 


195157 


11 


3 


196086 


196093 


2924373 


13 


3 


1682382 


1682389 


43820323 
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